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ABSTRACT 

We use a recently-developed analytic model for the ISM structure from scales of GMCs through 
star-forming cores to explore how the pre-stellar core mass function (CMF) and, by extrapolation, stellar 
initial mass function (IMF) should depend on both local and galactic properties. If the ISM is superson- 
ically turbulent, the statistical properties of the density field follow from the turbulent velocity spectrum, 
and the excursion set formalism can be applied to analytically calculate the mass function of collapsing 
cores on the smallest scales on which they are self-gravitating (non-fragmenting). Two parameters deter- 
mine the model: the disk-scale Mach number .M/, (which sets the shape of the CMF), and the absolute 
velocity /surface density (to assign an absolute scale). We show that, for normal variation in disk properties 
and gas temperatures in cores in the MW and local galaxies, there is almost no variation in the high-mass 
behavior of the CMF/IMF. The slope is always close to Salpeter down to < 1 Mq. We predict modest vari- 
ation in the sub-solar regime, mostly from variation in A4h, but this is consistent with the ~ lrr observed 
scatter in sub-solar IMFs in local regions. For fixed global (galaxy) properties, there is little variation in 
shape or "upper mass limit" with parent GMC mass or density. However, in extreme starbursts - ULIRGs 
and merging galaxy nuclei - we predict a bottom-heavy CMF. This agrees well with the IMF recently 
inferred for the centers of Virgo ellipticals, believed to have formed in such a process. The CMF is bottom 
heavy despite the gas temperature being an order of magnitude larger, because M.\ x is also much larger. 
Larger values make the "parent" cloud mass (the turbulent Jeans mass) larger, but promote fragmen- 
tation to smaller scales (set by the sonic mass, not the Jeans mass); this shifts the turnover mass and also 
steepens the slope of the low-mass CMF. The model may also predict a top-heavy CMF for the in-situ star 
formation in the sub-pc disk around Sgr A*, but the relevant input parameters are uncertain. 

Key words: star formation: general — galaxies: formation — galaxies: evolution — galaxies: active — 
cosmology: theory 



1 INTRODUCTION 

The origin of the stellar initial mass function (IMF) is a question 
of fundamental importance for the study of star formation, stellar 
evolution and feedback, and galaxy formation. It is an input into a 
huge range of models of all of these phenomena, and a necessary 
assumption when deriving physical parameters from many obser- 
vations. However, despite decades of theoretical study, it remains 
poorly understood. 

A particularly important question which has received great at- 
tention is whether the IMF is indeed "universal" or in fact varies in 
some environments. A growing body of observations find that the 
high-mass behavior of the IMF in the MW and nearby SMC/LMC 
systems is remarkably uniform, with a Salpeter-like (a = 2.3 in 
AN /AM oc M~ a ) slope consistently appearing in essentially all di- 
rectly measured systems: young clusters, globulars, disk, bulge, 
low/high mass GMCs, etc. (see e.g. |Kroupa|2002| |Chabrier|2003| 
|Bastia n et al. 2010 and references therein). There are sugges- 
tions of greater variation - albeit still modest - in how rapidly the 
IMF turns over in the sub-solar regime (< IMq), but the observa- 
tions are more uncertain at these masses (references above and e.g. 



A* may have a top-heavy IMF (either with a flat slope, a = 1.7 — 2, 
or a sharp sub-solar cutoff; see Paumard et al. 2006; Manes s~et al.| 
|2007|[Bartko et al.|2010) . These are believed to have formed in-situ, 
from a circum-BH thin disk (a quasi-accretion disk), as opposed to 
e.g. galactic center star clusters (that formed in their own clouds), 
and so may have had very different formation conditions from "nor- 
mal" stars. There have been a number of suggestions of top-heavy 
IMFs in high-redshift galaxies, from studies of galactic evolution, 
but these are all indirect, model-dependent constraints and subject 
to alternative interpretations (see Hopkins & Beacom 2006, van 
Dokkum|2008l|Dave|20"08l). Recently, though,|van Dokkum &Con- 



roy (2010 2011) have attempted to directly constrain the IMF in 



Chabri er|2005j|Allen et al.|2005| |Deacon et al.|2008| |Covey et al.| 
2008, De Marchi et al. 20101 



the centers of Virgo elliptical galaxies, and find the observations 
strongly favor a bottom-heavy IMF (a ~ 3). Since these galaxies 
formed their stars at high redshifts, this is a powerful constraint on 
any redshift evolution in the IMF. Moreover, it is well-established 
that stars in the central regions of ellipticals are formed via inflows 
of gas to large densities in extreme (probably merger-induced) nu- 
clear starbursts; this is the process seen "in action" in the centers 
of local UL IRGs (see|Kormendy & Sanders|1992||Mihos & Hern 
quist|1994||Hibbard & Yun|1999[ |Rothberg & Joseph|2004| |Ho~ 



There may be, however, variation in more extreme systems. 
The IMF in most of the galactic center is "normal" ( |Lockmann| 
|et al.|20 10 I, but there is compelling evidence that the pair of eccen- 
tric young circum-BH nuclear disks at just ~ 0. 1 — 0.5 pc from Sgr 
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kins ct al. [2008[ 1 2009| and references therein). So this should be 
directly related to the IMF in starburst environments, which many 
theoretical models have predicted should be top-heavy - the oppo- 
site of what is observed - owing to the higher observed gas tem- 
peratures and potentially large gas accretion rates (see e.g. Padoan 
|etal.|1997||Larson|2005) . 

Explaining the weak IMF variation in "normal" systems has 
actually been a challenge for theoretical models, which often pre- 
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diet large IMF variations with e.g. the local thermal Jeans mass 
(although see Jappsen et ah 2005 ). It is a particularly powerful con- 
straint for any model which attempts to predict IMF variations in 
systems where it has been suggested. There is no simple analytic 
model which has been shown to naturally explain both the weak 
IMF variation in MW and local group systems while also predict- 
ing large IMF variations in more extreme environments similar to 
those above. 

In this Letter, we therefore attempt to develop an analytic un- 
derstanding of IMF variation and apply it to these observations. It 
is increasingly clear that the ISM is governed by supersonic tur- 
bulence over a wide range of scales. Numerical and analytic argu- 
ments have shown that the density PDF in this regime approaches 
a lognormal distribution ( Vazquez-Semadeni 1994, Padoan ~et al.| 
[T9971 [QstrTker et al.|[l9"99| >. |Hennebelle & Chabrieij p008| > used 
this to propose a simple theoretical model for the stellar IMF, simi- 
lar to the derivation of the halo mass function in Press & Sc hechterl 
{1974}: if the PDF is lognormal, one can calculate the fraction of 
gas above a critical density (essentially the Jeans criterion for ther- 
mal+turbulent support) on a given scale, and associate this with the 
mass function. However, like the Press & Schechter ( 19741 argu- 
ment, this is only approximate as it does not resolve the "cloud in 
cloud" problem. Moreover, the focus therein was exclusively on 
small-scale properties, so a number of properties must be simply 
assumed rather than calculated from first principles and the shape 
and normalization of the IMF cannot be predicted a priori from 
global galaxy properties. Recently, Hopkins ( 201 lb) (hereafter Pa- 
per I) and Hopkins (2012) (Paper II) generalized this by showing 
how the mathematical excursion set formalism can be applied to 
turbulent density fields in galaxy disks, over the entire dynamic 
range from galactic scales to below the sonic length. This resolves 
these ambiguities, and in principle allows for a rigorous calculation 
of the core mass function (CMF), defined as the last-crossing distri- 
bution - specifically, the mass function of bound objects defined on 
the smallest scales on which they remain self-gravitating but do not 
have self-gravitating sub-regions (i.e. are not fragmenting). Paper 
II showed that the resulting CMF agrees well with CMFs observed, 
and (provided some plausible mean core-to-stellar mass conversion 
efficiency) can explain all of the major features of the canonical 
MW stellar IMF. Moreover, because the model self-consistently in- 
cludes the scales of the parent disk (where most of the turbulent 
power is concentrated), the distribution of parent cloud properties 
left as free parameters in Hennebelle & Chabrier ( 2008 ) can be de- 
rived; Paper I showed that the model is completely specified by the 
assumed turbulent velocity power spectrum and global disk prop- 
erties (surface density, velocity dispersion). 

In this Letter, we use the model for the core MF from Paper II 
to study how the predicted IMF might vary as a function of galactic 
properties and across different MW regions. In §[2] we summarize 
the model. In §[3] we derive some qualitative scalings for how the 
resulting CMF should vary with global parameters. In §[4] we com- 
pare the exact model results for different observed properties of 
the MW, starburst/ULIRG systems, and the circum-BH MW stel- 
lar disk. In §[5] we compare the model predictions for more typical 
variations in properties across the MW disk and in nearby galaxies, 
and examine how the CMF should vary within the MW in GMCs 
of different masses and regions of high/low densities. In § [6] we 
summarize our results and discuss the implications of the model. 

Before proceeding, we must caution that what the model here 
predicts directly is the pre-stellar core mass function, not the stellar 
IMF. It is by no means clear that the two can be trivially identi- 
fied with one another (modulo some mean "efficiency" of core-to- 



stellar mass conversion). In addition to the turbulent gas structure, 
the processes of fragmentation as cores contract and form proto- 
stars, accretion onto those protostars, outflows expelling material 
from cores, and ejection of protostars before accretion is exhausted, 
can all change the relation between CMF and IMF. It is beyond the 
scope of this paper to include all of these physics and therefore to 
make a rigorous absolute prediction of the IMF. Our intention is 
instead to study how the CMF varies as a function of global galaxy 
and turbulent properties. Our qualitative expectation is that other- 
wise identical CMFs, with similar microphysics, should produce a 
similar stellar IMF, and more bottom/top-heavy CMFs should cor- 
respond to more a bottom/top-heavy IMF. These assumptions are 
valid if the local star formation process depends only on e.g. the 
global core mass and/or stellar microphysics. But lacking a com- 
plete understanding of star formation, it remains possible that these 
microphysics will serve to "regulate" against variations in the IMF 
in a manner that offsets some of the predictions here. 

2 THE CMF & IMF 

If density fluctuations in supersonic turbulence are lognormal, then 
the variable 5(x) = ln[p(x) / po] + 5/2, where p(x) is the density at 
a point x, po is the global mean density and 5 is the variance in In p, 
is normally distributed according to the PDFQ 



./2+5 6XP ( 25, 

More generally, we can evaluate the field S(x\R), which is the 5(x) 
field averaged around the point x with some window function of 
characteristic radius R; this is also normally distributed, with a vari- 
ance S(R) that is directly related to the density power spectrum. 

Paper II showed that the CMF is equivalent to the last-crossing 
distribution, given by the numerical solution to the Volterra integral 
equation: 
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and 6(5) is the minimum value of the overdensity <5(x|6) which 
defines objects of interest (here, self-gravitating regions). 

In Paper I we derive S(R) and S(5) from simple theoretical 
considerations for all scales in a galactic disk. For a given turbulent 
power spectrum, with the assumption that the disk is stable (Q = 1), 
S(R) is determined by summing the contribution from the velocity 
variance on all scales R' > R 

3 vUk) 
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where W is the window function for the smoothing|^]6(i?) is 
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' The +S/2 term in S is required so that the integral of pPo(p) correctly 
gives po with (&) = 0. 

2 For convenience we take this to be a i-space tophat inside k < 1 /R, which 
is implicit in our previous derivation, but we show in Paper I and Paper II 
that this has little effect on our results. 
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where po is the mean midplane density of the disk, k = n/fl - 
for a constant-V,. disk, and 

a 2 g (R) = c 2 + {v 2 (R)) 



The mapping between radius and mass is 

M(*) S 4W[gr + (i + J)«p(-f 
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It is easy to see that on small scales, these scalings reduce to the 
Jeans criterion for a combination of thermal (c s ) and turbulent (v,) 
support, with M = (4ty/3) p CI i t R 3 ; on large scales it becomes the 
Toomre criterion with M — 7rE cr j t ^ 2 . The mass function is: 
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There are only two parameters that completely specify the 
model in dimensionless units. These are the spectral index p of 
the turbulent velocity spectrum, E(k) oc k~ p (usually p w 5/3 — 2), 
and its normalization, which we define by the Mach number on 
large scales Ml = (v}(h))/cj. The dimensional parameters h (or 
Cj) and po simply rescale the predictions to absolute units. 

In Hopkins (201 la I, we extended the derivation to solve the 
"two-barrier" problem. This is the solution for fe(S\5o[So]), i.e. 
given some value Sq on a larger-scale So < 5. We refer there for 
details, but note that this ultimately reduces to the solution for ft (S) 
above with the replacements B(S) -> B(S) - So, B(S') -> B(S') - 
5 ,S^S~- So, and 5' -»• 5" - So. 

3 QUALITATIVE SCALINGS 

For a choice of p and Mi,, it is straightforward to numerically 
determine the last-crossing mass function (CMF). Unfortunately a 
closed-form solution is not generally possible. However, we show 
in Paper II that on sufficiently small scales (near/below the sonic 
length), the "run" in S(R) ~ So(R SO mc) becomes small (since most 
of the power contributing in Eq.[5]comes from large scales), while 
B(S) rises rapidly, so dB/dS > B(S) /S > 1. In this limit, the MF 
can be approximated as 



p c (M) |dlnp c | 
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Since we are interested in the behavior on small scales, we 
drop higher-order terms in h and can re-write p cl -it in terms of the 
sonic length 



fl sonic = R[ V 2 , = c 2 ] = hM; 2/i "- l) 



which for p = 2 becomes just 
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with M w (4tv/3) p CI i t (R)R 3 . There is clearly a change in behav- 
ior below R « /? S onic- In Paper II we show that R > R m mc (where 
p cr it oc M _1//2 ) corresponds to the high-mass end of the CMF/IMF, 
giving a nearly power-law behavior, while R <C /? S onic (pcrit oc M~ ) 
corresponds to the low-mass end, where there is a quasi-lognormal 
turnover in the CMF and IMF. 

The turnover mass is related to M son i c = M(R S0T m), 



2 c'R., 
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where in the second equality we use Q ~ 1 (the standard Toomre Q 
parameter) to relate this to global parameters. 

Note that this is the global thermal Jeans mass in a disk - not, 
actually, the local Jeans mass oc G~ 3 ' 2 p -1 ' 2 , which is dimen- 
sionally the same but differs by powers of Mi, and p/ po. Both, 
however are different from the turbulent Jeans mass, 
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We showed in Paper I that M tlir b, j ea ns corresponds to the characteris- 
tic mass of GMCs. 

From Eq.[TT] when R 2> Z?sonic, the the high-mass slope of the 
CMF (dn/dM oc M~ a ) is approximately 
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where Mo = poh and the second term (in So) is small for the rea- 
sonable mass range (< 0.1), so we find 

3 
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(17) 



nearly independent of Mach number M h ■ The exact solution repro- 
duces this but with an even weaker p dependence (Paper II). 
On the other hand, the low-mass slope is 



aiowma SS « 3 - - (in 2 -4UM h - 2 In (M/M somc )) (18) 



which for typical parameters (using the approximate scaling So oc 
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4 CMF & IMF VARIATION IN EXTREME SYSTEMS 

We now consider an observationally motivated example. 

First, we consider canonical MW parameters: a minimum tem- 
perature T ~ 10K for the cores, disk gas surface density S gas ~ 
IOMq pc -2 , and large-scale turbulent velocity dispersion <r s {h) ~ 
10km s~' (Mi, ~ 30). For all cases, we will assume a spectrum 
p = 2, which is what is expected in highly super-sonic simulations 
and suggested by observations of the ISM (Burgers 1973; Larson] 
|1981| >. As discussed above and in Paper II, our results are not qual- 
itatively changed for p = 5/3 instead. Together this is sufficient to 
completely specify the model. 

Next, consider canonical parameters in ULIRG/starburst re- 
gions. Observations typically find T ~ 60 — 80 K in the dense 
molecular gas (n ~ 10 5 — 10 7 cm -3 ) with turbulent v elocity dis- 
persions a g (h) ~ 40- 100kms~' (Mi, ~ 50- 100) (|Downes & 



Solomon|19 98 ; Bryan t & Scoville|1999HWestmoquette et al.|2007 



Greve et al.|2009|>T and surface densities E gas ~ 10 — 10 Mq pc 
(e.g. |Kennicutt|1998) . 

In Fig. [T] we compare the CMF predicted in each case. Al- 
though the broad qualitative behavior is similar, there are striking 
differences. As expected from our analysis in §[3] the characteristic 
mass in the ULIRG case is somewhat smaller. From Eq. [14] com- 
pare M sor , ic ~ 1 .6M (MW) to Msonic ~ 0.4Mq (ULIRG). But even 
if we ignore the shift in the scale of M SO nic, the low-mass turnover 
occurs much more slowly in the ULIRG case, as expected from 
Eq.[T9]from the higher Mi,. The high-mass slope is nearly identical 
in both cases, as expected from Eq.JTT] but there is an intermediate- 
mass steepening in the ULIRG case that owes to second-order ef- 
fects from the run in S(R). 
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Figure 1. The exact predicted last-crossing mass function - i.e. star-forming 
core MF - from Eq.|2]for the properties of different observed systems. The 
predicted CMF mass scale is set by the sonic mass (Eq. |14) ; because of the 
assumption of marginal disk stability and a turbulent cascade, the shape is 
entirely specified by the dimensionless Mach number on disk scales Mi,. 
MW: Prediction for canonical local parameters: core minimum temperature 
T = 10K, disk gas surface density Eg as ~ IOMqpc -2 , and global tur- 
bulent velocity dispersion cr g {h) ~ lOkms -1 . We compare the observed 
|Kroupa|(2002) and [Chabner] {2003) IMFs, shifted in mass by an assumed 
star-to-core formation efficiency M* /M corc = 0.5. UL1RG: Prediction for 
properties typical in the central ~kpc of ULIRGs and mergers, the stellar 
remnants of which dominate the central light in most ellipticals: T ~ 65 K, 
S gas ~3x lO 3 M pc" 2 , OgO) ~ 80kms-'. We compare the bottom- 
heavy IMF (slope = —3 from 0.1 — 100Mq) fit to the centers of ellipticals 
in van Dokkum & Conroy 1 2010 1. BH: Prediction for properties of the MW 
~pc-scale circum-BH disk (these are less well-determined, see text). We 
compare a top-heavy IMF (slope = — 1.7; Bartko et al. 2010) over the pre- 
dicted range. 




-1 
Log[ M core /M e ] 

Figure 2. The CMF as Fig. [T] but for a more "normal" range of global 
parameters typical of the MW and local galaxies. We consider a few cases 
that vary crg as and £ gas together as expected for an exponential disk with 
V c =constant and Q = 1. We then consider independent variations of T, 
a, S gas . In all cases the high-mass behavior is nearly identical and close 
to Salpeter. The low-mass turnover still varies within this range, primarily 
with Mh, but less severely than in Fig. ^ We compare (grey shaded) the 
±1 a range of observed behavior from the IMF and PDMF compilations in 
|Bastian et ap(20T0) ; [De Marchi et ah](2l)T0l . 
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Figure 3. The CMF as Fig. ^ but for different "subregions" of the MW- 
like system. We compare the CMF in different isolated parent GMC masses 
(the parent largest mass/size scale on which the gas is self-gravitating, i.e. 
has (p(Scloud)) = Pcrit). for our standard parameters (solid lines; size and 
mass of the GMC labeled). For convenience, we normalize each model to 
the same maximum number density. We compare the IMF predicted in- 
side a parent cloud on the same spatial scale but with higher (dotted) or 
lower (dashed) mean density. Even though clouds have different cloud-scale 
Mach numbers ((.M 2 (i? c ioud)} increases as R\Q.a from ~ 0. I Mi, to 
~ 1.5.M/, in the smallest to largest clouds here) and density (p c rit/po de- 
creasing from ~ 50 to ~ 2.5), the CMF changes weakly with parent mass, 
provided it is large enough to properly sample the mass range. The high- 
mass end is suppressed in the lowest-mass clouds, but mostly because they 
lack sufficient mass to form many high-mass cores. 



An even more extreme example is the pc-scale circum-BH 
disk in the MW. Here, the initial gas conditions are uncertain; but 
from the observations, we can take the effective radius R ~ 0.5 pc 
and infer a surface density E gas ~ 0.6e - ' x 10 4 Mqpc -2 (where 
e = Mt/Mg as is the unknown star formation efficiency of the pro- 
genitor disk). Simulations of these circum-BH disks in Nayakshin 
et al. (2007); Hobbs & Nayakshin (2009) suggest an order-unity 
e ~ 0.5 and a % ~ 2km s~ , which is approximately what is needed 
to give Q— 1 (for Mbh = 4x 10 6 Mq, and k » f2 for the quasi- 
Keplerian orbit here). We adopt the same typical nuclear gas tem- 
peratures as the ULIRG case. The CMF in this case is very dif- 
ferent from the previous cases. This is because the system is ex- 
tremely thin and « large, dominated by the external BH potential, 
so Mi, is relatively small, but also our assumptions that /fgonic "C h 
and dB/dS 2> 1 break down and the the correction terms in ft not 
captured in Eq.|l l|dominate. 

5 CMF & IMF VARIATION IN DIFFERENT LOCAL 
REGIONS 

The above examples are intentionally extreme. Normal galaxies ex- 
hibit a much more narrow range of E gas , T, and a g . In Fig. [2] we 
compare the CMF predicted for a "typical" range in these parame- 
ters. In all these cases, the high-mass slope is nearly identical and 
close to the Salpeter value. Modest changes in Mi, still lead to dif- 
ferences at the low-mass end, although the shape is quite similar. 
We compare these variations to the observed range inferred for the 
IMF at these masses in different systems (or equivalently, at these 
masses, the PDMF); the range in "turnover speed" is consistent. 

There is also the question of how the CMF/IMF should vary 
within different GMCs or sub-regions of the Galaxy with fixed 
global properties. To address this we use the solution to the two- 
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barrier problem, which lets us calculate the CMF for any given 
overdensity on a larger "parent" scale. In Paper I, we identified self- 
gravitating regions on larger scales (which contain sub-regions of 
fragmentation), i.e. crossing the barrier B(Sq) at smaller 5o <S(M), 
as GMCs. We can therefore calculate the CMF for a GMC of 
size/mass scale Ro (So) or Mo (So) by taking the two-barrier solution 
with So(So) = B(Sq). We show this in Fig. [5]for a range of "par- 
ent" GMC scales. We can further arbitrarily vary So on each scale: 
we consider regions more/less dense than the "collapse threshold" 
(typical GMC) at the same spatial scale. The results are typically 
less sensitive to parent cloud properties than the other parameters 
we have considered. Because the global parameters are fixed, par- 
ent scale does not change Mi, or the sonic length - it is largely 
equivalent to shifting the "background" mass normalization. In the 
lowest-mass GMCs, the high-mass end is suppressed, but this oc- 
curs only when the GMC masses are so low that they simply cannot 
contain the highest-mass cores. 

6 DISCUSSION 

We have used the method from Paper I and Paper II, which devel- 
oped a formal means to calculate the statistical properties of bound 
objects formed by turbulent density fluctuations in a galactic disk, 
to study the predicted variation in the pre-stellar core MF, and by 
implication the stellar IMF, across different environments. Specif- 
ically, in Paper I, we showed how the full excursion set formalism 
can be applied to lognormal density distributions in a galactic disk, 
and demonstrated that, for a given turbulent power spectrum shape 
E(k) oc k~ p , their dimensionless properties were only a function 
of the normalization of that power spectrum on large scales M (h). 
In Paper II, we extended this to define the "last crossing distribu- 
tion" which defines the mass spectrum of self-gravitating cores on 
the smallest scales on which they are bound; following Hennebelle 
|& Cha brier ( 2008 I we showed that for typical conditions, this pre- 
dicts a mass function in good agreement with observed CMFs and 
(for plausible core-to-stellar mass conversion efficiencies) consis- 
tent with the shape of the canonical stellar IMF. 

Here, we show that under "normal" conditions - specifically, 
a range of star-forming core temperatures T ~ 5 — 20 K, galaxy 
velocity dispersions a g ~ 5 — 15kms~', and galaxy surface den- 
sities S gas ~ 1 — 20Mq pc -2 - there is almost no variation in the 
high-mass (Salpeter) behavior. In Paper I we showed this is also 
true if we vary the turbulent power spectrum shape p = 5/3 — 2. 
The low-mass end is somewhat sensitive to the global Mach num- 
ber Mi,, turning over more slowly in systems with higher Mi, (per 
Eq. [T9| >. The same result was obtained with a more approximate 
derivation of the CMF in |Hennebelle & Chabrier| ( |2008l[2009l ). The 
change is not especially large, however, corresponding to a range of 
slopes a ~ — 1 at ~ 0.1 Mq (with larger, but much more uncer- 
tain, variation below this). We stress that this is the global Mh, it 
does not necessarily imply variation in regions with different local 
Mi, (since their density field has contributions from all scales). 

To address that question, in Hopkins (2011a), we extend the 
excursion set model to derive the solution to the two-barrier last- 
crossing problem. This corresponds to the solution for the CMF in 
sub-regions of any specified size/mass scale, within a common par- 
ent system of fixed global properties. We apply this here to examine 
how the CMF should vary in different GMC sub-regions within the 
MW. Since in Paper I we show that GMCs can be identified as so- 
lutions to the first-crossing problem, we can specifically show how 
the CMF should vary in "typical" GMCs of different masses, and 
higher/lower density regions. The variation is minimal in all regions 



of modest mass (despite systematically different densities and local 
M). 

At low GMC masses the CMF is truncated, but this is basi- 
cally just because there is not enough collapsing mass to form the 
most massive cores. Since the mass density in GMCs is primarily in 
high-mass systems and this occurs only in the lowest-mass clouds, 
it will have negligible global effects. The mass scale of the pre- 
dicted CMF does not vary systematically with parent GMC mass. 

These predictions agree with the fact that there is remarkably 
little variation in the high-mass IMF behavior observed in different 
MW and local group star forming regions. And the predicted low- 
mass variation in these regions in fact quite similar to the scatter 
inferred between different regions observationally (although some 
of this may owe to observational uncertainties). 

If, however, we consider much more extreme systems, we pre- 
dict significant variations in the CMF. For the typical conditions in 
ULIRGs and other dense, starburst regions, we predict a bottom- 
heavy CMF. There is strong observational and theoretical evidence 
that the stars in the central ~kpc of essentially all ellipticals must 
have formed in nuclear (probably merger-induced) starbursts, and 
the properties of recent "ULIRG/starburst relics" agree well with 
elliptical/bulge centers (see e.g. Kormendy & Sanders 1992; Hib- 
bard & Yun|1999|[Rothberg & Joseph 2004| |Hopkins et al.|2008| 



2009). In fact Hopkins & Hernquist (2010) showed that the cen- 



tral densities and mass profile structure of ellipticals imply, if they 
obeyed any sort of Kennicutt-Schmidt relation in their formation, 
that they must have been starbursts with densities and velocity dis- 
persions quite similar to what we assume in Fig.[T] This may, there- 
fore, explain recent observational suggestions that the central re- 
gions of nearby Virgo ellipticals exhibit a bottom-heavy IMF simi- 
lar to that predicted here ( |van Dokkum & Conroy|2010|[20TT) . 

The key reason for this is that these regions have much larger 
global Mach numbers Mi, than the typical MW regions, even if we 
allow for much larger minimum gas temperatures T ~ 50 — 100 K. 
Mathematically, the predicted bottom-heavy CMF in ULIRGs cor- 
responds to the slow low-mass turnover predicted in Eq. [19] al- 
though there are higher-order corrections that only appear in the 
exact solution (not the approximate CMF derivation in Hennebelle 
|& Cha brier ( 2008)). Physically, this is a result of turbulent fragmen- 
tation. Although the characteristic turbulent Jeans mass - which de- 
termines the largest structures - increases (scaling oc a 4 ), this ap- 
pears in the ^r.s?-crossing distribution: i.e. it means bigger GMCs, 
not bigger cores. The GMC/first-crossing structures are not directly 
relevant, since they will be self-gravitating on many smaller scales 
- i.e. fragment - internally. But the CMF turnover is not exactly 
set by the local thermal Jeans mass either. Strictly, the mass scale 
of the last-crossing distribution is set by the sonic mass, so there 
is some competition between higher Mi, (which drives the sonic 
length down) and higher c s . If the system is globally stable, we can 
equate this mass scale to purely global parameters, and infer that 
the rise in surface density E gas is sufficiently rapid that it "wins 
out" and lowers the sonic mass. However, even if this mass scale 
were the same or somewhat higher, the CMF would still be quite 
bottom-heavy: a higher Mh not only shifts the sonic mass scale, 
but also imprints larger density fluctuations in the turbulent field 
on all scales. This promotes more fragmentation over a more broad 
range of masses, "slowing down" the low-mass turnover or even 
steepening the slope towards low masses. 

This is an example where the distinction between the first and 
last-crossing distributions is critical: if we ignored it, we would ar- 
rive at similar conclusions to Pad oan et al.| l|l997 1 - essentially, the 
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exact opposite of what we find - where the authors argued that star- 
burst systems should have more top-heavy IMFs. 

Another extreme system is the ~ 0.5pc-scale circum-BH disk 
in the MW, which observations have suggested may have a top- 
heavy IMF < |Paumard et al.|2006[|Maness et al.|2007||Bartko et al.| 
|2010| >. We attempt to extend our estimate to this system as well, al- 
though we caution that the characteristic parameters during its for- 
mation are quite uncertain and our model may not be applicable at 
all because of non-linear tidal and eccentricity corrections ( |Alexan-| 
der et al. 2008 ). The predicted CMF is somewhat top-heavy, but it 
also has an unusual shape which is much more restricted in mass 
range (truncated by shear and tidal forces). Part of this top-heavy 
character owes to its having a relatively small Mi, if it has 2=1 
(for a Q = 1 disk near a dominant BH, a g /V c ~ M^/Mbh ~ 10~~ 
here). But it is complicated because the disk has such a small scale 
height that the distinction between first and last crossings is blurred. 
Nevertheless, it is suggestive and agrees with results from numer- 
ical hydrodynamic simulations (Nayakshin et al. 2007; H obbs &| 
|Nayakshin|2"009) . 

We caution that, as discussed in Hennebelle & Chabrier 
(20091, the low-mass end is also the transsonic regime (below the 
sonic length), so our simple assumption of isothermal gas will 
break down at some level. This is probably the largest uncertainty 
in our model; a complicated equation-of-state, magnetic fields, and 
radiative feedback are certainly capable of modifying the relation 
between CMF and IMF, especially if we were to account for time- 
dependent fragmentation in contracting cores, and our predictions 
should not be taken too literally as an absolute prediction for the 
IMF (especially at low masses). Rather, our intention is to moti- 
vate some potential observed IMF variations (or lack thereof) as a 
consequence of the properties of the pre-stellar core mass function 
with global galaxy properties. 
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